version February 25, 2013: fm 

The WISE gamma-ray strip parametrization: 
the nature of the gamma-ray Active Galactic Nuclei of Uncertain 

type 

(N 

F. Massaro 1 , R. D'Abrusco 2 , G. Tosti 3 - 4 , M. Ajello 1 , D. Gasparrini 5 , J. E. Grindlay 2 & 

^vq Howard A. Smith 2 . 

$_^ 

SLAC National Laboratory and Kavli Institute for Particle Astrophysics and Cosmology, 

2575 Sand Hill Road, Menlo Park, CA 94025 
Harvard - Smithsonian Astrophysical Observatory, 60 Garden Street, Cambridge, MA 02138 



\o 



P_] Dipartimento di Fisica, Universita degli Studi di Perugia, 06123 Perugia, Italy 

X 

Istituto Nazionale di Fisica Nucleare, Sezione di Perugia, 06123 Perugia, Italy 



£ 



ASI Science Data Center, ESRIN, 1-00044 Frascati, Italy 



43 
Oh 

6 

ABSTRACT 

03 

Despite the large number of discoveries made recently by Fermi, the origins 
of the so called unidentified 7-ray sources remain unknown. The large number 



> 

O 

CO of these sources suggests that among them there could be a population that sig- 

nificantly contributes to the isotropic gamma-ray background and is therefore 

^ crucial to understand their nature. The first step toward a complete comprehen- 

sion of the unidentified 7-ray source population is to identify those that can be 
associated with blazars, the most numerous class of extragalactic sources in the 

> 7-ray sky. Recently, we discovered that blazars can be recognized and separated 

from other extragalactic sources using the infrared (IR) WISE satellite colors. 
The blazar population delineates a remarkable and distinctive region of the IR 
color-color space, the WISE blazar strip. In particular, the subregion delineated 
by the 7-ray emitting blazars is even narrower and we named it as the WISE 
Gamma-ray Strip (WGS). In this paper we parametrize the WGS on the basis of 
a single parameter s that we then use to determine if 7-ray Active Galactic Nuclei 
of the uncertain type (AGUs) detected by Fermi are consistent with the WGS 
and so can be considered blazar candidates. We find that 54 AGUs out of a set 
60 analyzed have IR colors consistent with the WGS; only 6 AGUs are outliers. 
This result implies that a very high percentage (i.e., in this sample about 90%) 
of the AGUs detected by Fermi are indeed blazar candidates. 



Subject headings: galaxies: active - galaxies: BL Lacertae objects - radiation 
mechanisms: non-thermal 



1. Introduction 

With the recent advent of the Fermi mission, the 7-ray astronomy is living a new golden 
age with several striking discoveries already performed during the first three years. 

According to the second Fermi 7-ray LAT catalog (2FGL, lAbdo et al. 20111) . Fermi 
detected 1873 sources, 576 of which are still unidentified even if the localization of the 7-ray 
sources has significantly improved with respect to the past 7-ray missions. For this reason, 
despite the large number of new discoveries already achieved, the nature of the unidentified 
7-ray sources is still an open question. This unsolved issue is extremely relevant for the origin 
of the isotropic gamma-ray background, since, given the large number of unidentified 7-ray 
sources, new classes of unknown extragalactic 7-ray sources that can significantly contribute 
to the isotropic gamma-ray background could be hidden. 

On the other hand, the most detected 7-ray sources in the MeV-GeV energy range be- 
long to the rarest class of active galactic nuclei, the blazars. They are an intriguing class of 
active galactic nuclei, characterized by non-thermal radiation emitted over the entire electro- 
magnetic spectrum, and interpreted as arising from a relativistic jet closely aligned to the line 
of sight (see e.g. IBlandford fc Rees 1978|) . Blazars come in two flavors: the BL Lac objects 



and the flat spectrum radio quasars, where the common discriminating criterion between the 
two classes is the equivalent width of the optical emission lines, traditionally weaker than 
5 A in the former rather than in the latter (IStickel et al. 199"T| IStoke et al. 199 ill . In the 



following, we indicate the BL Lacs as BZBs and the Flat Spectrum radio quasars as BZQs, 
according to the ROMA-BZCAT nomenclature (IMassaro et al. 2009| IMassaro et al. 2010J) . 



Recently, using the preliminary data release of the WISE infrared (IR) survey (Wright et al. 2010) 
we discovered that IR color-color diagrams allows us to distinguish between extragalac- 
tic sources dominated by non-thermal emission, like blazars, and other classes of galaxies 
and/or active galactic nuclei (Massaro et al. 2011, hereinafter Paper I, see also Plotkin et al. 
2011). In particular, the blazar population delineates a tight, distinct region of the IR color 
space, indicated as the WISE Blazar Strip (Paper I). The WISE Blazar Strip is a region 
in the 3D infrared color space delineated by the blazar population. This region is narrower 
when considering only the IR colors of the blazars that are detected in the 7-rays, indicated 
as the WISE gamma-ray strip (WGS, see ID'Abrusco et al. 20121 hereinafter Paper II). A 
3D scatter plot of the WISE Blazar Strip and the subregion of the WGS are shown in the 



IR diagram of Figure [II while the [3.4]-[4.6]-[12] /zm 2D projection is reported in Figure [2J 

One of the major difficulties of the association procedures for the Fermi 7-ray sources 
with active galactic nuclei is that, due to the lack of radio and X-ray informations and to 
the large uncertainty on the 7-ray position, it is not always possible to recognize if there is 
a blazar candidate within the positional error region. Thus, the main aim of this paper is 
to build a parametrization of the WGS in order to verify whenever 7-ray sources have been 
associated to a counterpart that is a blazar candidate, being consistent with the WGS. In 
particular, we studied the WISE counterparts of the Active galactic nuclei of uncertain type 
(AGUs), defined according to the 2FGL and the 2LAC criteria (lAckermann et al. 201 ip . 



and their consistency with the WGS. The AGUs are defined as the radio and/or X-ray 
counterparts of 7-ray sources associated to by the 2FGL Likelihood Ratio method, but 
without a good optical spectrum that enable their classification (lAbdo et al. 201 ill . This 
paper is organized as follows. Section [2] describes the procedure adopted to parametrize the 
WGS. Section [3] we discuss on the consistency of the WISE counterparts of the AGU sample 
with the WGSwhile Section [4] describes the non-parametric analysis of the WGS based on 
the Kernel Density Estimation (KDE). In Section [5] we investigated possible selection effects 
that could affect our WGS parametrization. Our conclusions are discussed in Section |6l 



2. The WISE gamma-ray strip parametrization 

2.1. The sample selection 

We use the sample of 7-ray emitting blazars already selected in Paper II to parametrize 
the WGS. 

This sample was selected from the 2FGL, that contains 805 sources associated with a 
blazar: 435 BZBs and 370 BZQs respectively. However only 659 (347 BZBs and 312 BZQs) 
of these are listed and classified according to the criteria used in the ROMA-BZCAT (e.g., 
IMassaro et al. 20091) . We excluded from our analysis all the blazars with a Fermi analysis 



flag, according to the 2FGL and the 2LAC flAbdo et al. 201 H lAckermann et al. 20111 ). 

In particular, 329 (164 BZBs and 165 BZQs) blazars, out of the original 659, lie in 
the portion of the sky reported in the WISE preliminary source catalog, but only 296 (143 
BZBs and 153 BZQs) have a WISE counterpart within 2.4" radius (see Paper I). To be more 
conservative, we excluded from our analysis 12 blazars (8 BZBs and 4 BZQs) with respect to 
the 296 blazars in sample selected in Paper II, because they have a 95% upper limit on the 
WISE magnitude at 22 /zm. Then, we use this 2FB sample composed of the 284 blazars (135 
BZBs and 149 BZQs) to build the WGS parametrization. We notice that all the selected 



blazars also belong the 2LAC sample (lAckermann et al. 201 ip . 



According to the classification available in ROMA-BZCAT the blazars of uncertain type 
have been excluded from our analysis, while the BL Lac candidates have been considered as 
BZBs. More details about the 2FB sample and the source selections are given in Paper II. 

Finally, we emphasize that our selection is based only on 7-ray blazars that belong to 
the ROMA-BZCAT because this is the largest catalog of blazars available in literature in 
which each source is spectroscopically classified at optical frequencies. 



2.2. The WISE blazar associations 

The IR color-color diagrams have been built using the archival WISE Preliminary Source 
Catalog, that covers ~ 57% of the sky Q . The WISE mission mapped the sky at 3.4, 4.6, 
12, and 22 fim in 2010 with an angular resolution of 6.1, 6.4, 6.5 & 12.0" in the four bands, 
achieving 5a point source sensitivities of 0.08, 0.11, 1 and 6 mJy respectively in unconfused 
regions on the ecliptic. All the WISE magnitudes are in the Vega system. In particular, the 
absolute (radial) differences between WISE source-peaks and "true" astrometric positions 
anywhere on the sky are no larger than ~ 0.50, 0.26, 0.26, and 1.4" in the four WISE bands, 
respectively ( ICutri et al. 201 1q 



For our analysis, unless stated otherwise, we considered only WISE sources detected 
with a minimum signal-to-noise ratio of 7 in at least one band[j The positional coincidences 
of blazars in the observed WISE sky have been searched within a circular region of radius 
2.4". This corresponds to the combination of the error of 1", assumed for the radio position 
reported in the ROMA-BZCAT (IMassaro et al. 2"0 09) and taking into account of astrometric 



uncertainties in the WISE preliminary data, and the positional error of the fourth WISE 
band at 22/xm (i.e., 1.4") (see also Paper I). All the associations of the 2FB blazars with 
WISE sources are unique and no multiple matches have been found (see Papers I and II 
for more details). The chance probabilities of the WISE associations for the sources in the 
2FGL and in the ROMA-BZCAT are reported in Paper II. 



1 wise2.ipac.caltech.edu/docs/release/prelim/preview. html 

2 wise2. ipac.caltech.edu/docs/release/prelim/expsup/sec2_3g. html 

3 We take the opportunity to correct here an error that appears in Paper II. The sources of the 2FB 
sample are detected with a minimum signal-to-noise ratio of 7 in at least one band rather than in all four 
bands as reported in Paper II. 



2.3. The WISE gamma-ray strip projections 

We built the parametrization of the WGS considering only the sources in the 2FB sample 
(see Section 2.1[) and using the three different 2D projections of the WGS delineated in the 



[3.4]-[4.6]-[12] /im, [4.6]-[12]-[22] /mi, [3.4]-[4.6]-[22] /xm color-color planes. In each color-color 
2D projection we determined the smaller irregular quadrilateral containing at least 95% of 



the blazars in the 2FB sample considering their position within error (see also Section 2.4 
for more details). The irregular quadrilateral defining the WGS subregions have been drawn 
by hand. The KDE analysis has been used to verify, a posteriori, that the hand-drawn 
boundaries of the WGS are in agreement with the sharp decline in density of WGS sources, 
as evaluated by this non parametric method (see Section |4J for more details). This WGS 
modeling has been developed separately for the BZBs and the BZQs, and for all their three 
2D projections. 

The [3.4]-[4.6]-[12] /im, [4.6]-[12]-[22] /im and [3.4]-[4.6]-[12]-[22] /im 2D projections of 
the WGS for the BZB and the BZQ populations are shown in Figure |3j Figure |4j Figure [5] 
and Figure [6]), respectively. 

In the following we also report the boundaries chosen for our WGS parametrization. 
For the BZB projections the extremal points of the WGS have coordinates: Pj=(2. 01,0.37), 
P 2 =(3.30,1.17), P 3 =(2.59,1.20), P 4 =(1.52,0.51) in Figure[3) Pi = (2.20,1.65), P 2 =(2.72,2.57), 
P 3 =(2.29,3.30), P 4 =(1.20,1.96), in Figure @ (upper panel), while Pi = (2.05,0.33), P 2 =(2.83,1.07), 
P 3 =(2.28,1.21), P 4 =(l. 20,0.73), in Figure [B] (lower panel). 



On the other hand, for the BZQ projections the extremal points of the WGS have 
coordinates: Pi = (2.90,0.85), P 2 =(3.81,1.17), P 3 =(3.29,1.67), P 4 =(2.29,1.08) in Figure gj 
Pi = (2.25,2.22), P 2 =(3.04,3.05), P 3 =(2.67,3.70), P 4 =(l. 68,2.85), in Figure § (upper panel), 
while P 1 =(2.48,0.78), P 2 =(3.05,1.17), P 3 =(2.55,1.50), P 4 =(1.72,1.12), in Figure § (lower 
panel) . 



2.4. The strip parameter s 

To illustrate the WGS parametrization we consider the schematic case of the first pro- 
jection: [3.4]-[4.6], [4.6]-[12], hereinafter Ci-c 2 with the correspondent errors <7i and <r 2 , re- 
spectively (see Figure [7|). 



Based on the WISE source location in the C1-C2 diagram we can distinguish 5 types of 
objects. Each source, given its IR colors, corresponds to a single point in each 2D color-color 
projection of the WGS. However, including the errors on both axes, it is represented by a 



cross with 4 extremal points, calculated considering the ±lcr error on each color. Then, we 
can define five different type of sources, according to the schematic view shown in Figure [7| 

• type 4'- sources with all the extremal points within the WGS projection; 

• type 3: sources for which only 3 extremal points lie within the region of the WGS; 

• type 2: sources with only two extremal points consistent with the WGS; 

• type 1: sources with only a single extremal point associated with the WGS; 

• type 0: sources without extremal points of the error cross on the WGS. 

We can assign to each type of source a discrete strip parameter d ranging between and 
1, according to the scheme illustrated in Figure [7} For example, in the case of the c\-c 2 
projection of the WGS, we assign to type 4 sources a value d 12 =l, while source of type 
corresponds to d 12 =0. For the same 2D projection, the intermediate values have been 
assigned as follows: type 3 have <ii2=0.75, type 2 have <ii2=0.5 and type 1 di 2 =0.25. 

On the same c\-c 2 diagram, we also assign a weight strip parameter wu to each value 
of the d 12 parameter defined as: W\ 2 = (o"i 02) -1 , proportional to the area of the ellipse 
described by the error bars of each point (see inset of Figure [7] for more details). Then, we 
define the continuos strip parameter s\ 2 as: 

s 12 = d 12 w 12 . (1) 

We note that the parameter 1012 has been chosen to take into account the different errors on 
both axes when comparing two sources that might belong to the same type. It also allow us 
to make the Si 2 continuos rather than discrete as d 12 . 

We repeated the entire procedure described above for each 2D projection of the WGS: 
C1-C2, C2-C3 and C1-C3, generating the values of the continuos strip parameters si 2 , S23, S13, 
respectively. Then, all these values of the strip parameters for three different 2D projections 
have been combined together to define an unique total strip parameter s. The total strip 
parameter is the geometric average of the s values of each 2D projection: 

s = (s 12 s 23 s 13 ) 1/3 . (2) 

We emphasize that sources that lie outside of the WGS in at least one of its 2D projection 
have one of the correspondent si 2 , s 23 , S13 parameter equal to zero and consequently the total 
s value is null as well. This occurs because the discrete d parameter is zero for sources outside 
the WGS (see Figure [7}. 



We divided all the s parameters for the maximum s values of the BZBs and BZQs that 
lie on the WGS to re-normalize s in the range and 1. This re-normalization can be applied 
to the s values of all the WISE sources, because those outside the WGS will have s null by 
definition. 

The s parameter represents an estimate of the distance, in the IR colors parameter space 
and weighted with the errors on each axes, between the WGS and a generic WISE source, 
that could potentially belong to it; this s parameter is different from zero only in the case 
in which the all color error bars of a WISE source are consistent with the WGS. Therefore 
these s values can be used to rank each IR WISE source according to their association to 
the WGS. 

Finally, we note that to test if a generic WISE source has IR colors consistent with the 
BZBs or with the BZQs subregion of the WGS, the total strip parameters are indicated as 
Sb and s q , respectively. We introduced the above divisions for the s parameters because in 
future works this allow us to verify if a generic WISE source is more consistent to be a BZB 
or a BZQ, enabling a classification for new IR sources that could lie on the WGSbeing 7-ray 
blazar candidates. 



2.5. The s parameter distributions 

We considered a sample composed of all the WISE sources lying in two circular regions 
of 1 deg radius, centered at high and low Galactic latitude b, with the center coordinates of 
(7, 6) = (255,-55) deg and (7, 6) = (338,-l) deg, respectively. These sources do not have upper 
limits on their WISE magnitude values and are detected with a signal to noise ratio > 7 in 
at least one band as for the blazars in the 2FB sample. 

We calculated the s parameters for all the 11599 generic IR WISE sources. This analysis 
provides an estimate of the probability to find a generic WISE source in the sky with a 
particular value of Sb and/or s q . We perform this analysis considering the distinction between 
the two blazar classes (i.e., BZBs and BZQs), The distributions of the Sb and s q parameters 
for the BZBs and the BZQs that lie on the WGS in comparison with the generic IR WISE 
sources are shown in Figure [8] and in Figure |9j respectively. 

From the distributions of the Sb and s q parameters for the generic IR WISE sources, 
we note that 99.9% of them have s b <0.24 and s q <0.38. Then, for the BZBs in the 2FB 
sample only 6 sources out of 135 have Sb < 0.24, and in the case of the BZQs only 33 sources 
out of 149 show s q values lower than 0.38. We also note that 99.0% of the generic IR WISE 
sources have Sb <0.10 and only 2 BZBs are below this value, while 97.2% of the generic IR 



WISE sources together with only 5 BZQs out of 149 have s q <0.14. Finally, on the basis of 
the above s distributions we define the outliers of the WGS, WISE sources that have values 
of s b <0.10 or s q <0.14. 

We recognize that the above choice of the s& and s q thresholds are extremely conser- 
vative. This choice has been made on the basis of the actual sky coverage of the WISE 
preliminary data release. At the present status of our analysis we are not able to investigate 
the IR emission of all the blazars that are listed in the 2FGL, and the 2FB sample used to 
build the WGS parametrization is small with respect to how it will be available when the 
WISE full archive will be released. Consequently, regarding the choice of the threshold values 
for the s parameters, we preferred the efficiency to the completeness of our method selecting 
the s limiting values from their distributions in low galactic latitude regions, even if this 
choice could increase the possible contamination of the WGS. A deeper investigation of this 
problem will be considered in future as an a posteriori analysis of WGS parametrization. In 
particular, when the WISE full release will be available, we will improve our method taking 
into account of the IR source density at different galactic latitudes and of the varying depth 
of the exposure for the WISE observations. 



3. The AGU counterparts on the WGS 

We considered the sample composed of all the AGUs already classified in the analysis 
of the 2FGL and the 2LAC flAbdo et al. 20111 lAckermann et al. 20111 ). According to the 
2FGL, the AGUs could be all blazar candidates without a good optical spectrum or without 
an optical spectrum at all ( lAbdo et al. 20lT| lAckermann et al. 20111) . 

We selected the AGUs that lie in the portion of the sky surveyed by WISE during the 
first year, corresponding to 148 Fermi sources. Then, we excluded from our analysis all the 
AGUs with a Fermi analysis flag, according to the 2FGL and the 2LAC. 

The association between each AGU counterpart and the WISE sources have been evalu- 
ated on the basis of the same criterion chosen for the blazars on the WISE Blazar Strip (see 



Section 2.2 and Paper I for more details), considering the position of the radio counterpart 
for each AGU as reported in the 2FGL and/or in the 2LAC. There are 60 AGUs out of 148 
for which there is a unique association with a WISE source (see Section [3]) within the usual 
region of 2.4" radius and with a chance probability of 0.008, estimated adopting the method 
described in Maselli et al. (2010) (see also IMaselli et al. 20l"H and Paper I) and without 
upper limits on the WISE magnitudes within the WISE preliminary data release. 

Subsequently, we used the IR colors of the AGU counterparts, as associated in the 



2LAC, to verify if the WISE counterparts of the 7-ray sources in the 60 AGU sample lie on 
the WGS, evaluating their s values following the procedure described in Section [2} 

We found that 6 outliers do not belong to the WGS out of 60 AGUs, according to 
the inclusion based on the threshold values of s& <0.10 or s q <0.14. With this analysis on 
the 60 AGU sample, we have been able to check if the association provided by the 2FGL 
corresponds to a blazar lying on the WGS. 

Finally, we estimated the IR spectral index am using the [3.4]- [4.6] /zm color according 
to Eq. (1) of Paper II, and we evaluated the correlation between am and the spectral index 
of the associated 2FGL source a 7 . We found a linear correlation between am and ct 7 for the 
54 AGU that lie on the WGS, with a correlation coefficient p= 0.56 and a chance probability 
of 8.96 x!0~ 6 and a slope m=0.30±0.06, that is consistent with that of the WGS blazars 



(p=0.68, m=0. 36±0. 02, see Paper II) within one sigma (see Figure 10). On the other hand, 
the 6 outliers have a weaker linear correlation between the two spectral indices then the 
previous sample with p=0.40 (chance probability of 0.08) and m=0.12±0.07, different from 
that of the blazars on the WGS (see Figure 10). Finally, in Table [I] and in Table [2] we report 
the colors, the IR spectral indices and the s parameters together with the 2FGL name and 
the WISE and the counterpart names of each AGU analyzed. The class of each AGU as 



derived from the 2LAC analysis is also indicated (lAckermann et al. 20111) . 



4. An independent non-parametric analysis: the kernel density estimation 

To test our analysis, we also performed a statistical investigation based on an inde- 
pendent non-parametric method as the KDE technique as already proposed in Paper I (see 
also ID'Abrusco et al. 2009( ILaurino &: D'Abrusco 20TT| and reference therein). The KDE 
method provides an effective way of estimating the probability function of a multivariate 
distribution and do not require any assumption about the shape of the "parent" distri- 



butions. In Figure 11, the isodensity contours drawn from the KDE density probabilities 
and associated with different levels of density are plotted for the blazars of the WGS in its 

[3.4]-[4.6]-[12]-[22] /im 2D projection. 

Consequently, for a generic source in the WISE archive we can provide an estimate of 
the probability 7Tfcrf e that a blazar of the WGS has the same IR colors, this is a surrogate of 



the probability that a WISE source is consistent with the WGS. In Figure [TTJ we also show 
the AGU counterparts with respect to the isodensity contours of the WGS, to highlight the 
outliers. We also report in Table [T] and in Table [2j the value of nkde f° r each AGU analyzed. 

Finally, we note that there are some AGUs for which the KDE analysis suggests that 
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the source is not consistent with the WGS, even if the parametric method indicates it as a 
possible candidate. The reason for this to happen is that, as previously mentioned, the KDE 
method does not take into account the errors on the IR colors. As a consequence, sources far 
from the WGS but with large errors could be associated to low density values, as calculated 
by the KDE method, and discarded. However, our parametrization of the WGS allows us 
to take into account the errors on the WISE colors and to consider also this type of sources. 
We emphasize that all the sources that the WGS parametrization indicates as outliers have 
also nkde typically lower than ~ 1% of being consistent with the WGS. 



5. An analysis of possible selection effects 

In future, thanks to the developed WGS parametrization, we will be also able to inves- 
tigate if there are selection effects that could affect our analysis as for example driving the 
WGS thickness. At the current stage of our study, we are able to estimate when a generic 
WISE sources in consistent with the WGS itself that is a necessary tool to compare different 
samples for future investigations; Thanks to the parametrization developed we will be also 
able to verify if there are IR blazars that belong to the WGS but are not detected in 7-rays, 
and which could be the physical conditions if this occurs. 

We remark that the link between the IR and the 7-ray properties of blazars, is mainly 
due to the relation between the blazar spectral shape in the IR and in the 7-rays. To evaluate 
if selection effects due to flux limits in the selected sample could affect our WGS description 
we performed the following tests. We restricted our analysis to the bright WISE blazars with 
IR magnitudes in the ranges: mi <13.5, m 2 <12, m 3 <11, m 4 <7.5, that belong to both the 
WISE Blazar strip and to the WGS, and we found that a difference in their thickness is still 



evident (see Figure 12 for the standard 2D projection in [3.4]- [4.6]- [12] /im color diagram). 



This plot suggests that the origin of the WGS is not due to a selection of bright IR blazars. 

We also compared the WGS as formed by the blazars present in the 2FGL and those 
detected in the first Fermi LAT catalog (1FGL) (lAbdo et al. 20101) . and again we did not 
find any clear difference between the WGS drawn with the bright or the faint 7-ray blazars. 
This again suggests that the relation between the WGS and the 7-ray detectability is related 



to the blazar spectral shape (see Figure 13 for the standard 2D projection in [3.4]- [4.6]- [12] 
^m color diagram). On the basis of the WGS parametrization presented in this paper, the 
above issues will be deeply addressed in future works, after the full release of the WISE 
all-sky survey. 



11 



6. Summary and Discussion 

On the basis of the recent results on the characterization of the IR colors of blazars 
provided by WISE (Paper I) and on the comparison with their 7-ray emission (Paper II), we 
developed a method based the WISE Blazar Strip to identify blazar counterparts of 7-ray 
sources. 

We developed a method to parametrize the WGS in the 3D color diagrams based on 
its 2D projections. This method is characterized by the use a continuos parameter s, in the 
range 0-1, that takes into account of the errors of all the IR colors and provides clues 
on the position of a generic WISE source relative to the WGS in the 3D color parameter 
space. High values of the s parameters are associated to sources that lie inside the WGS (as 
the 7-ray blazar population of the 2FB sample). The WGS has been parametrized in two 
subregions, the first containing the BZBs while the other with the BZQs although in the 
present work we are only interest in searching for blazar counterparts that lie on the WGS. 

We applied our parametrization to the sample of the AGUs selected from the 2FGL and 
the 2LAC. We found that there are 148 AGUs that can be analyzed within the footprint of 
the WISE preliminary source catalog. However, according to our association procedure (see 



Secton 2.2) only 60 AGUs have a unique WISE counterpart without any upper limit on the 



WISE magnitude values. 

Then, we calculated the distributions of their s parameter and found that 54 out of 
60 AGUs analyzed are consistent with the WGS, corresponding to the 90% of the 7-ray 
counterparts analyzed while the remaining 6 AGU counterparts are outliers of the WGS. 
In particular, for the 54 AGUs that are consistent with the WGS we also found that the 
correlation between the otm and a 1 is in agreement with that found for the blazars that 
constitute the WGS itself while the same correlation for the 6 outliers is inconsistent with 
it. 

We also applied the KDE non-parametric test to obtain the probability that an AGU 
counterpart belong to the WGS and we found consistent results with our parametrization 
(see Section 14] for more details). 



In addition, an extensive investigation of all the unidentified 7-ray sources in the 2FGL 
that fall in the area of the sky where the WISE preliminary data have been already released 
will be provided in a forthcoming paper ( IMassaro et al. 2012ap . Searching for blazar candi- 
dates within the unidentified 7-ray source sample could potentially leading to the discovery 
new class of 7-ray emitting sources. 

Further improvements of the WGS parametrization will be also possible in the future, 
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when the whole WISE catalog will be available and this parametric method would be cal- 
ibrated at different b values not only to look for counterparts of 7-ray sources but also to 
search for new blazar candidates all over the sky (IMassaro et al. 2012b}) . 
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Fig. 1. — The 3D representation of the WISE Blazar Strip (blazars are indicated in magenta) 
and the subregion of the WGS (7-ray emitting blazars are indicated in black) in the IR colors 
built with the magnitudes in the WISE bands at [3.4]-[4.6]-[12]-[22] fan. 



a 




[4.6] - [12] (mag) 



Fig. 2. — The 2D projection of the WISE Blazar Strip (blazars are in magenta) and the 
subregion of the WGS (7-ray emitting blazars are in black) in the IR color diagram [3.4]- 
[4.6]-[12] /im. 
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BZB population: [3.4]-[4.6]-[12] projection 




2.1 2.4 2.7 3 3.3 3.6 3.9 

[4.6] - [12] 



Fig. 3. — The [3.4]-[4.6]-[12] /zm 2D projection of the WGS in the subregion of the BZB 
population is shown. 



BZQ population: [3.4]-[4.6]-[12] projection 




[4.6] - [12] (mag) 



Fig. 4. — Same of Figure [| but for the case of the BZQ population. 



BZB population: [4 .6] -[12] -[22] projection 



BZB population: [3.4]-[4.6]-[12]-[22] projection 
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Fig. 5. — Same of Figure [3] for the case of the BZB population on the WGS in the two 
remaining different color-color projections. 



BZQ population: [4.6]-[12]-[22] projection 



BZQ population: [3.4]-[4.6]-[12]-[22] projection 
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Fig. 6. — Same of Figure [4] for the case of the BZQ population on the WGS in the two 
remaining different color-color projections. 
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The strip parameter s 12 association 
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Fig. 7. — The schematic view of the strip parametrization in the example of the c\-Ci 2D 



projection in arbitrary units (a.u.). We report the method described in Section 2.4 to assign 



to each point of the 2D projection of the WGS a value of the discrete strip parameter d 12 
and the associated value of the weight strip parameter W\i- The combination of these two 
values provides the continuos strip parameter si 2 for each given source (see Equation [T]). 
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Fig. 8. — The distribution of the strip parameter Sb for the BZBs (blue) and the BZQs (red) 
that lie on the WGS in comparison with the generic IR WISE sources (black). 
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Fig. 9. — Same of Figure [8| for the distribution of the s q parameter. 




Fig. 10. — The correlation between am and a 7 for the WISE counterparts of the AGUs. 
The background grey circles represent the correlation found for the WGS blazars (see Paper 
II), while the red filled squares are the AGUs that have been found consistent with the WGS 
accordingly with our parametrization. The remaining 6 outliers (blue filled circles) show a 
weaker and different correlation than the other samples. 
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Fig. 11. — The isodensity contours drawn from the KDE technique for the blazars of the 
WGS (grey circles) are shown for the case of the [3.4]-[4.6]-[12]-[22] /zm 2D projection. The 
AGU identified as blazar candidates (filled black squares) are also shown in comparison with 
the outliers (black open square) to show the consistency between the WGS parametrization 
and the KDE analysis. 
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Fig. 12. — The 2D projection of the WISE Blazar Strip (magenta) and the subregion of the 
WGS (black) in the IR color diagram [3.4]-[4.6]-[12] /xm when only bright IR blazars are 
considered (i.e., those with WISE magnitudes in the ranges: mi <13.5, m 2 <12, m 3 <11, 
1714 <7.5, respctively). 
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Fig. 13. — The 2D projection of the WGS in the IR color diagram [3.4]-[4.6]-[12] /im when 
only bright 7-ray blazars listed in the 1FGL (black) are considered in comparison with those 
in the 2FGL (red). There are no clear differences on the thickness of the WGS in this 2D 
projection due to the different samples considered. 



